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How much helicity is needed to drive large-scale dynamos? 
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Magnetic field generation on scales large compared with the scale of the turbulent eddies is known to be 
possible via the so-called a effect when the turbulence is helical and if the domain is large enough for the 
a effect to dominate over turbulent diffusion. Using three-dimensional turbulence simulations, we show that 
the energy of the resulting mean magnetic field of the saturated state increases linearly with the product of 
normalized helicity and the ratio of domain scale to eddy scale, provided this product exceeds a critical value 
around unity. This implies that large-scale dynamo action commences when the normalized helicity is larger 
than the inverse scale ratio. Recent findings of a smaller minimal helicity may be due to the onset of small- 
scale dynamo action at large magnetic Reynolds numbers. The onset of large-scale dynamo action is difficult to 
establish when the kinetic helicity is small. 

PACS numbers: 



I. INTRODUCTION 

The origin of magnetic fields in astrophysical bodies like 
the Earth, the Sun and galaxies is studied in the field of dy- 
namo theory. The temporal variation and strength of those 
fields rules out a primordial origin, through which the mag- 
netic field would have been created in the early Universe. For 
magnetic fields with energies of the equipartition value, i.e. 
the turbulent kinetic energy of the medium, the primordial hy- 
pothesis explains their strength after creation, but falls short of 
explaining how the field outlives billions of years of resistive 
decay QJJ] . 

In dynamo theory, astrophysical plasmas are considered 
sufficiently well conducting fluids where the inertia of the 
charge-carrying particles can be neglected. In this approxi- 
mation the equations of magnetohydrodynamics (MHD) pro- 
vide an adequate model of the medium. In this framework it 
has been studied under which conditions magnetic fields of 
equipartition strength and scales larger than the turbulent mo- 
tions are created and sustained |01 . 

A successful theoretical model describing the dynamo's be- 
havior is the mean-field theory. It relates the small-scale tur- 
bulent motions to the mean magnetic field via the so-called a 
effect, which provides the energy input via helical turbulent 
forcing. During the kinematic phase, i.e. negligible back re- 
action of the magnetic field on the fluid, the a effect gives a 
positive feedback on the large-scale magnetic field, which re- 
sults in its exponential growth. Already the consideration of 
the kinematic MHD equations with negligible Lorentz force 
sheds light on the growth rate of the different modes of the 
magnetic field during the kinematic phase. In the kinematic 
phase the growth rate A at wave number fc is given by [2] 



A = ak — rfrk 2 = (C a — 1)^t^ 2 , 



(1) 



where C a — a/rjTk is the relevant dynamo number for the 
a 2 dynamo, a is the a coefficient which is proportional to 
the small-scale kinetic helicity, and ryx = V + % is the sum 
of molecular and turbulent magnetic diffusivity. Clearly, dy- 
namo action occurs when \C a \ > C™\ where the onset 
condition is C£ rlt = 1, Standard estimates for isotropic 



turbulence in the high conductivity limit fH [U yield a = 
— (t/3)(u-u) and?7t = (t/3)(m 2 ), where r is the correlation 
time of the turbulence, u> = V x u is the vorticity and u is 
the velocity in the small-scale fields. Here (.) denote volume 
averages. Using r/ t rj, we have 
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(w • u)/k(u 2 



(2) 



It is convenient to define (u-u)/ k{ (it 2 ) as the normalized ki- 
netic helicity, ef, so C a — —eiki/k. This scaling implies that 
the critical value of the normalized helicity ef scales inversely 
proportional to the scale separation ratio, i.e. t{ oc (fcf/fc) -1 , 
where k -C fcf is the wave number of the resulting large- 
scale magnetic field. This wave number can be equal to 
k = k\ = 2-k/L, which is the smallest wave number in a 
periodic domain of size L. 

In summary, the critical dynamo number C„ rlt , which de- 
cides between growing or decaying solutions of the large- 
scale dynamo (LSD), is proportional to the product of nor- 
malized helicity ef and scale separation ratio fcf jk. Therefore, 
the amount of helicity needed for the LSD is inversely propor- 
tional to the scale separation ratio, and not some higher power 
of it. It should be noted that the normalized kinetic helicity 
ef used here is not the same as the relative kinetic helicity, 
if = (u> ■ u) / (w rms u rms ). The two are related to each other 
via the relation 



ef/ef = ku/kt, 



(3) 



where k w ~ ui Tms /u Ims is inversely proportional to the Tay- 
lor microscale. Here, the subscripts rms refer to root-mean- 
square values. For small Reynolds numbers, fc w provides a 
useful estimate of the wave number fcf of the energy-carrying 
eddies. In contrast, for large Reynolds numbers Re, we expect 
k u /k[ to be proportional to Re 1 / 2 , so e~f decreases correspond- 
ingly while ef remains unchanged. Below we shall present 
new results suggesting that in the present case of helical tur- 
bulence the exponent is not 1/2 but 1/4. 

To understand the saturation of a helical dynamo, it is im- 
portant to understand the relation between the resulting large- 
scale field and the associated small-scale field. Indeed, the 
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growth of the large-scale field is always accompanied by a 
growth of small-scale magnetic field. Small-scale here means 
the scale of the underlying turbulent motions, which drive 
the dynamo. Conservation of total magnetic helicity causes 
a build up of magnetic helicity at large scales and of oppo- 
site sign at small scales [4, 5]. As the dynamo saturates, the 
largest scales of the magnetic field become even larger, which 
finally leads to a field of a scale that is similar to that of the 
system itself. This can be understood as being the result of 
an inverse cascade, which was first predicted based on closure 
calculations |6]. 

If the domain is closed or periodic, the build up of small- 
scale magnetic helicity causes the a effect to diminish, which 
marks the end of the exponential growth and could occur well 
before final saturation is reached. The dynamo then is said 
to be catastrophically quenched and, in a closed or periodic 
system, the subsequent growth to the final state happens not 
on a dynamical timescale, but on a resistive one. Quenching 
becomes stronger as the magnetic Reynolds number increases, 
which, for astrophysically relevant problems, means a total 
loss of the LSD within the timescales of interest. In the case 
of open boundaries magnetic helicity fluxes can occur, which 
can alleviate the quenching and allow for fast saturation of the 
large-scale magnetic field |7|-ll0ll. 

In a recent publication ifTlll it was argued that for periodic 
boundaries the critical value of ef for LSD action to occur de- 
creases with the scale separation ratio like eji rlt oc (fcf//ci)~ 3 . 
Their finding, however, is at variance with the predictions 
made using equation ([T), which would rather suggest a depen- 
dence of e"' 1 * oc (fcf/fci) -1 with C™* = 1. This discrepancy 
could be a consequence of the criterion used in [11] for de- 
termining C" 1 *. The authors looked at the growth rate of the 
magnetic field after the end of the kinematic growth phase, but 
only at a small fraction of the resistive time. Therefore their 
results might well be contaminated by magnetic fields result- 
ing from the small-scale dynamo (SSD). Earlier simulations 
IU2I1 have demonstrated that for ReM > 100, the growth rate 
of the helical LSD approaches the well-known scaling of the 
non-helical SSD with A oc Re 1 / 2 , which corres pon ds to the 
turnover rate of the smallest turbulent eddies 

Given that the LSD is best seen in the nonlinear regime 
lfl5tl . we decided to determine C„ rlt from a bifurcation dia- 
gram by extrapolating to zero. In a bifurcation diagram we 
plot the energy of the mean or large-scale field versus C a . 
Simple considerations using the magnetic helicity equation 
applied to a homogeneous system in the steady state show 
that the current helicity must vanish lfl5ll . In a helically driven 
system, this implies that the current helicity of the large-scale 
field must then be equal to minus the current helicity of the 
small-scale field. For a helical magnetic field this yields the 

value of (B } / B 2 q approximately as etkt/ki, i.e., propor- 
tional to C a . Here B cq = (^op) 1 ^ 2 Wrms is the equipartition 
value of the magnetic field with the permeability /in. and the 
mean density p. Again, this suggests that the LSD is excited 
for ef > (fcf/fci) -1 rather than some higher power of k(/ki. 
This is a basic prediction that has been obtained from nonlin- 
ear mean-field dynamo models that incorporate magnetic he- 
licity evolution [16] as well as from direct numerical simula- 



tions in the presence of shear lfl7ll . The purpose of the present 
paper is to reinvestigate these results over a broader parameter 
range in the light of recent conflicting findings Ultl . 

II. THE MODEL 

A. Basic equations 

The underlying equations for non-relativistic MHD, which 
we solve numerically, are given as 

^-A = U x B - Wo J, (4) 
at 

^-U = -c 2 Vlnp+-JxB + F visc + f, (5) 
Dt p 

£lnp=-V-C/ ) (6) 

where A is the magnetic vector potential, U the velocity, B 
the magnetic field, 77 the molecular magnetic diffusivity, po 
the vacuum permeability, J the electric current density, c s 
the isothermal sound speed, p the density, -F v isc the viscous 
force, / the helical forcing term and D/D< = d/dt + U ■ V 
the advective time derivative. The viscous force is given as 
-Fvisc = p~ x V • 2vpS, with the traceless rate of strain tensor 
S with components Sy = + — ■ U and the 

kinematic viscosity v. Commas denote partial derivatives. 

The energy supply for a helically driven dynamo is pro- 
vided by the forcing function / = f(x, t), which is a helical 
function that is random in time. It is defined as 

/(as, t) = Rc{Nf k{t) exp[ifc(i) • x + ty(t)]}, (7) 

where x is the position vector. The wave vector k(t) and the 
random phase — tt < <fi(t) < tt change at every time step, 
so f(x,t) is ^-correlated in time. For the time-integrated 
forcing function to be independent of the length of the time 
step St, the normalization factor N has to be proportional 
to 5t~ x / 2 . On dimensional grounds it is chosen to be N = 
/oCsd/clcg/^) 1 / 2 , where /o is a non-dimensional forcing am- 
plitude. We choose fo — 0.02 which results in a maximum 
Mach number of about 0.3 and an rms value of about 0.085. 
At each timestep we select randomly one of many possible 
wave vectors in a certain range around a given forcing wave 
number. The average wave number is referred to as kf . Trans- 
verse helical waves are produced via II 1411 

fe = R • f k ' with Rtj = J ■, (8) 

V 1 + c 

where a is a measure of the helicity of the forcing and a = 1 
for positive maximum helicity of the forcing function. Fur- 
thermore, 

/<» ohcl > = (fe x e) /v/fc 2 - (fe • e) 2 (9) 

is a non-helical forcing function, where e is an arbitrary unit 
vector not aligned with k; note that \fk\ 2 = 1 and 

/fe-(ifex f k )* =2ak/(l + a 2 ), (10) 
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so the relative helicity of the forcing function in real space is 
2*/(l + a 2 ). 

For comparison with earlier work, we shall also use in one 
case an ABC-flow forcing function, 



Jo 



sin X% + a cos Xi 
sin X\+a cos X 3 



v/|(l + ct 2 ) ysin X 2 + a cos X Y 



(11) 



where Xi = kfXi + 6% and 6i = 9o cos Uit are time-dependent 
phases that vary sinusoidally with frequencies caj and ampli- 
tude So- This forcing function is easy to implement and serves 
therefore as a proxy of helical turbulence; see Ref. EMU, 
where the phases changed randomly. We have restricted our- 
selves to the special case where the coefficients in front of 
the trigonometric functions are unity, but those could be made 
time-dependent too; see Ref. [19]. However, as we will see 
below, ABC-flow driven dynamos do not show some crucial 
aspects of random plane wave-forced helical turbulence. Most 
of the results presented below concern the forcing function (O 
and only one case with Eq. (fTTT > will be considered at the end. 

Our model is governed by several non-dimensional parame- 
ters. In addition to the scale separation ratio kf/ki, introduced 
above, there are the magnetic Reynolds and Prandtl numbers 



Re 



M 



(12) 



These two numbers also define the fluid Reynolds number, 
Re = u xms /i/k{ = Rcm/Pim- The maximum values that 
can be attained are limited by the numerical resolution and 
become more restrictive at larger scale separation. The cal- 
culations have been performed using the PENCIL Code (see 
http://pencil-code.googlecode.com) at resolutions of up to 
512 3 mesh points. 



B. Mean-field interpretation 

The induced small-scale motions u are helical and give rise 
to the kinetic a effect yfl 



a K = 



[u> ■ u) 

3w rms fcf 



(13) 



The relevant a effect for dynamo action is the sum of the ki- 
netic and magnetic a 1I20I1 : 



u) + (j-b)/(p) 
3u rms fcf 



(14) 



The resulting large-scale fields are partially helical, which 
means one can write 



(J-B)=e m k m (B z 



(15) 



with large-scale wave vector k m and corresponding fractional 
helicity e m . However, in the cases considered below the do- 
main is triply periodic, so the solutions are Beltrami fields for 
which k m — ki and e m = 1 is an excellent approximation, 
and only ef will take values less unity. Nevertheless, in some 



expressions we retain the e m factor for clarity. For example, 
the saturation value of the large-scale magnetic field, B sat , is 
given by Hill 



B Lt/ B L - (\C a \/e„ 



Ik 



(16) 



where C a = ock/ {ijrki) is the relevant dynamo number 
based on the smallest wavenumber in the domain and 1 
I + 3/ReM = ifi/rit is a correction factor resulting from the 
fact that r/T is slightly bigger than rjt. The factor 3 in the ex- 
pression for l results from our definition of ReM and the fact 
thatJIl 

Vt = Wrms/3fcf = 7?Re M /3. (17) 

Equation ( [ToT l shows clearly the onset condition \C a \ > 
|e m | fa 1. Using Eqs. O and (17), we find 



(u> • u) 



ik\ 



(18) 



ik\ M 2 ms 

From equation ( [ToT l we can derive the critical value of the nor- 
malized helicity ef as a function of the scale separation ratio 
Setting C a to is critical value (\C a \ = e,„) we obtain 



(19) 



which is at variance with the findings in 111 111 . 

Once the dynamo is excited and has reached a steady state, 
not only a but also ?/ t will be suppressed. This can be taken 
into account using a quenching factor g(B), so ijt(B) = 
VtQfj(B) with g = (1 + g|B|/_B cq ) _1 . Equation ( fTBT ) is then 
modified and reads S 2 at /_B 2 q = (\C a \ — C a o)i/e m with 

C a0 = [1 - (1 - g)/i]e m . (20) 

Note that C a o = ej,, 1 in the unquenched case, i.e., for g = 1, 

C. Simulation strategy 

We recall that our forcing term / in equation (01 is a 
stochastic forcing centered around the wave number k{. In 
contrast to 111 111 this forcing is delta-correlated in time. The 
fractional helicity of the helical forcing is a free parameter. 
The simulation domain is a periodic cube with dimensions 2ir. 
Due to the cubic geometry of the domain the large-scale mag- 
netic field can orient itself in 3 possible directions. Therefore 
we compute three possible planar averages (xy, xz, and yz 
averages). From their resistive evolution we infer their satu- 
ration values at the end of the resistive phase. The strongest 
field gives then the relevant mean-field B. 

Since B is helical and magnetic helicity can only change 
on resistive timescales, the temporal evolution of the energy 
of the mean magnetic field, M(t), is given by lfl5Tl 



M(t) = M Q - M 1 e~ t/T , 



(21) 



-B 2 at is the square 



where r _1 = 27]ef n kf is known, M 
of the desired saturation field strength, and Mi is an un- 
known constant that can be positive or negative, depending 
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FIG. 1: Example showing the evolution of the normalized (B ) 

2 2 

(dashed) and that of {B ) + rd{B ) /At (dotted), compared with its 
average in the interval 1.2 < 2i)k\t < 3.5 (horizontal blue solid 
line), as well as averages over 3 subintervals (horizontal red dashed 
lines). Here, B is evaluated as an xz average, (B) xz . For compari- 
son we also show the other two averages, {B) xy (solid) and {B) xy 
(dash-dotted), but their values are very small. 
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on whether the initial magnetic field of a given calculation 
was smaller or larger than the final value. (Here, an initial 
field could refer to the last snapshot of another calculation 
with similar parameters, for example.) The functional behav- 
ior given by Eq. $2l[ allows us to determine B 2 &t as the time 
average of M + rdM/dt, which should only fluctuate about 
a constant value, i.e., 



B 



1 



t 2 -h 



dt'. (22) 



This technique has the advantage that we do not need to 
wait until the field reaches its final saturation field strength. 
Error bars can be estimated by computing this average for 
each third of the full time series and taking the largest depar- 
ture from the average over the full time series. An example 

is shown in Fig. Q] where we see (B ) still growing while 

2 2 2 

(B )+rd(B ) /dt is nearly constant when (B ) has reached 
a value less than half its final one. 



III. RESULTS 



A. Dependence of kinetic helicity on a 



FIG. 2: Dependence of relative kinetic helicity ei (a) and normal- 
ized kinetic helicity ef (b) on the helicity parameter a of the forcing 
function Eq. $Q together with the analytical expression 2<r/(l + a 2 ) 
(solid line). 
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We recall that the relative helicity of the forcing function is 
/ • V x //[/ rms (V x /) rms ] = 2a/(l + a 2 ). This imposes 
then a similar variation onto the relative kinetic helicity, if — 
(u>-u)/w rms u rms ; see Fig.|2|a). However, as discussed above, 
if is smaller than ef by a factor k u /kf, which depends on the 
Reynolds number. It turns out that ef matches almost exactly 
the values of 2a/(l + a 2 ); see Fig.Efb). 

The theoretically expected scaling k u /kf cx Re 1 / 2 is a 
well-known result for non-helical turbulence 0231, but for 



FIG. 3: Dependence of k^/kt on Re. The open and closed circles 
correspond to runs with Ptm = 1, while the + signs correspond to 
runs with Ptm = 100, so Re = Reni/PrM. 



helical turbulence those scaling relations are modified; see, 
e.g., Ref. j23ll . For our current data we find that k u /kf ~ 
0.8 Re 1 / 4 ; see Fig. [3] The reason for this power law behavior 
is however unclear. 
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FIG. 4: Steady state values of (B ) / as a function of C a to- 
gether with the theoretical prediction from equation j 1 6b (dashed 
line) and a linear fit (dotted line). 



B. Dependence on scale separation 

We perform simulations with different forcing wave num- 
bers kf and different values of ef at approximately constant 
magnetic Reynolds number, ReM ~ 6, and fixed magnetic 
Prandtl number, PrM = 1. Near the end of the resistive sat- 
uration phase we look at the energy of the strongest mode at 
k = ki, using the method described in § 111 CI 

As seen from Eq. ( [ToT l. mean-field considerations predict 
a linear increase of the saturation magnetic energy with C a 
and onset at C a = 1. This behavior is reproduced in our 
simulation (Fig. |4j), where we compare the theoretical predic- 
tion with the simulation results. For different values of kf/k\ 
and C a we extrapolate the critical value C„ rlt « 1.2 (Fig.|4]i, 
which gives the critical values ej; rlt s» 1.2i(fcf/fci) = 
1.7 (fcf/fci) -1 for which the LSD is excited. For each scale 

separation value we plot the dependence of (B ) /B^ q on ef 
(Fig. [5]l and make linear fits. From these fits we can extrap- 
olate the critical values e™*, for which the LSD gets excited 
(Fig. [6}, which gives again e™* sa 1.7 (fef/fci) -1 . 

It is noteworthy that the graph of (B ) / £? 2 q versus C a de- 
viates systematically (although only by a small amount) from 
the theoretically expected value, [C a — l)c. While the slope is 
rather close to the expected one, the onset of LSD is slightly 
delayed and occurs at C a « 1.2 instead of 1. The reason for 
this is not clear, although one might speculate that it could be 
modelled by adopting modified effective values of i or e m in 
Eq. ( f20b . Apart from such minor discrepancies with respect to 
the simple theory, the agreement is quite remarkable. Never- 
theless, we must ask ourselves whether this agreement persists 
for larger values of the magnetic Reynolds number. This will 
be addressed in § |III CI 

At this point we should note that there is also a theoretical 
prediction for the energy in the magnetic fluctuations, namely 
(b 2 )/Bg « (Cq-C^VCq. Nonetheless, the results shown 



FIG. 5: Steady state values of (B }/B^ q as a function of ef for 
various scale separation values fef/fei together with linear fits. 
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FIG. 6: Critical value for the normalized kinetic helicity ef for which 
LSD action occurs for different scale separations. 

in Fig.|7]deviate from this relation and are better described by 
a modified formula 

(6 2 )/5 e 2 q oc 1 - (Cf 7<7 Q )" (withn«4). (23) 

Again, the reason for this departure is currently unclear. 

C. Dependence on Re M 

To examine whether there is any unexpected dependence of 
the onset and the energy of the mean magnetic field on ReM 
and to approach the parameters used in 111 ill , who used values 
up to ReM = 1500, we now consider larger values of the mag- 
netic Reynolds number. This widens the inertial range signifi- 
cantly and leads to the excitation of the SSD. Higher values of 
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FIG. 7: Steady state values of (b 2 ) /B 2 q as a function of Ca together 
with the fit formula from equation i23i with n = 4, compared with 
n = 1 (dotted) and n = 2 (dashed). 
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FIG. 8: Steady state values of (B ) / B 2 ^ as a function of C a for 
PrM = 100 (dashed line) and PrM = 1 (solid line) for k{ /ki = 5 
and different values of ReM (different symbols), compared with the 
theoretical prediction (dotted line). 



ReM can more easily be reached at larger values of PrM- This 
is because at large values of PrM, most of the injected energy 
is dissipated viscously rather than resistively, leaving less en- 
ergy to be channelled down the magnetic cascade [24]. This is 
similar to the case of small values of PrM, where larger fluid 
Reynolds numbers can be reached because then most of the 
energy is dissipated resistively [12]]. Here, however, we shall 
be concerned with the former case of large values of PrM- 

The results for PrM = 100 are presented in Fig. [8]for differ- 
ent values of ef and ReM being either 80, 200, or 600. Most 
importantly, it turns out that the critical value for the onset of 
LSD is not much changed. We now find C™* « 0.9 instead 




FIG. 9: Visualization of B x on the periphery of the domain for 
PrM = 100 after resistive saturation. 

of 1. Furthermore, the dependence of (B ) / B^ q on C Q , and 
thus also the value of C£ rlt , are the same for all three values 
of ReM- It is however surprising that it is now even below 
unity. Furthermore, the slope is different from that of Fig. [4] 
it is now closer to 2 than to the earlier value of 1 . This dis- 
crepancy with the theory can be easily explained by arguing 
that the relevant value of B cq has been underestimated in the 
low PrM cases. Looking at the power spectrum of the high 
PrM simulations in Fig. [TUf a), we see that the kinetic energy 
is indeed subdominant and does not provide a good estimate 
of the magnetic energy of the small-scale field (6 2 )/2/io- By 
contrast, for Ptm = 1, the magnetic and kinetic energy spec- 
tra are similar at all scales except near k = k\; see Fig. fTOf b). 
The slight super-equipartition for k > kf is also typical of 
small-scale dynamos [ 14]. 

A visualization of the magnetic field in this case is given 
in Fig.|£l where we show B x on the periphery of the compu- 
tational domain. The magnetic field has now clearly strong 
gradients locally, while still being otherwise dominated by a 
large-scale component at k = k\. In this case the large-scale 
field shows variations only in the y direction and is of the form 

B = (sin kiy, 0, cos k\y) B sat . (24) 

This field has negative magnetic helicity, so J ■ B = —kiB 2 , 
as expected for a forcing function with negative helicity. 
We have argued that the reason for the larger slope in the 

2 

graph of (B ) versus C a is related to B cq being underesti- 
mated for large values of Pi - m- To confirm this, we now con- 
sider calculations with PrM = 1, two different values of ef, 
and fixed values of ReM (either 160 or 320), and fixed scale 
separation ratio kf/ki = 5. We see in Fig.[8]that the slope is 
indeed smaller. Furthermore, C™' is now above 1 , and even 
larger than at low ReM (now C™ 1 * sa 1.5 instead of 1.2). 
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FIG. 10: Comparison of kinetic and magnetic energy spectra for 
PrM = 100 (upper panel) and PrM = 1 (lower panel) for a — 0.2 
(solid lines) and 0.12 (dashed lines). 




2 4 6 
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FIG. 11: Similar to Fig.Q] but for time-dependent ABC-flow driv- 
ing. As in Fig.Q] we have here fcf/fci = 15 and ReM ~ 4. 



D. ABC-flow forcing 

In this paper we have used the fact that the saturation field 
strength is described by Eq. ( fTol l. While this is indeed well 
obeyed for our randomly driven flows, this does not seem to 
be the case for turbulence driven by ABC-flow forcing. We 
demonstrate this by considering a case that is similar to that 
shown in Fig. [TJ where ReM ~ 6 in the saturated state. We 
thus use Eq. ( fTTT ) with a = 8 = 1 and fcf/fci = 15. The 
kinematic flow velocity reaches an equilibrium rms velocity 
of Uq — fo/vkf. The magnetic Reynolds number based on 



this velocity is Uo/r]kf, which is chosen to be 13, so that dur- 
ing saturation the resulting value of ReM is about 6, just as in 
Fig.Q] For the x, y, and z components we take different forc- 
ing frequencies such that coj/fcii/o is 10, 11, and 9 for i = 1, 
2, and 3, respectively. These values correspond approximately 
to the inverse correlation times used in Hill . The result is 
shown in Fig. QT| It turns out that the magnetic field grows 
initially as expected, based on Eq. d2Tb . but then the final sat- 
uration phase is cut short below B^ at /B^ q 3 rather than the 
value 12 found with random wave forcing. This is reminis- 
cent of inhomogeneous dynamos in which magnetic helicity 
fluxes operate. In homogeneous systems, however, magnetic 
helicity flux divergences have only been seen if there is also 
shear B25I1 . In any case, the present behavior is unexpected and 
suggests that the effective value of C a is reduced. Using the 
test-field method [26, 27], we have confirmed that the actual 
value of C a is not reduced. The dynamo is therefore excited, 
but the value implied for the effective helicity is reduced. 

Another possibility is that, especially for small values of a, 
the ABC-flow has non-generic dynamo properties that emu- 
late aspects of large-scale dynamos. An example is shown in 
Fig.[T2]where we plot the time evolution of all three planar av- 
erages (yz, xz, and xy). Even for a — 0.01, large-scale mag- 
netic fields are still excited, but the field orientation changes 
periodically on a timescale of 1-2 diffusion times. This is 
obviously a fascinating topic for further research, but it is un- 
related to our main question regarding the minimal helicity of 
generic turbulent dynamos. The main point of this section is 
to emphasize the limited usefulness of ABC-flow dynamos. 
Another such example are dynamos driven by the Galloway- 
Proctor flow, which also has a number of peculiar features; see 
Ref. JH. 
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IV. CONCLUSIONS 

In this paper we have studied the simplest possible LSD 
and have investigated the dependence of its saturation ampli- 
tude on the amount of kinetic helicity in the system. We recall 
that the case of a periodic domain has already been investi- 
gated in some detail [29, 30], and that theoretical predictions 
in the case with shear [16] have been verified numerically for 
fractional helicities ltl7ll . Yet the issue has now attracted new 
interest in view of recent results suggesting that, in the limit of 
large scale separation, the amount of kinetic helicity needed to 
drive the LSD might actually be much smaller than what ear- 
lier calculations have suggested. This was surprising given 
the earlier confirmations of the theory. As explained above, 
the reason for the conflicting earlier results may be the fact 
that the LSD cannot be safely isolated in the linear regime, 
because it will be dominated by the SSD or, in the case of the 
ABC-flow dynamo, by some other kind of dynamo that is not 
due to the a effect. While there will always remain some un- 
certainty regarding the application to the much more extreme 
astrophysical parameter regime, we can now rule out the pos- 
sibility of surprising effects within certain limits of ReM be- 
low 600 and scale separation ratios below 80. 

In stars and galaxies, the scale separation ratio is difficult 
to estimate, but it is hardly above the largest value considered 
here. This ratio is largest in the top layers of the solar con- 
vection zone where the correlation length of the turbulence is 
short (1 Mm) compared with the spatial extent of the system 



(100 Mm). 

Of course, the magnetic Reynolds numbers in the Sun and 
in galaxies are much larger than what will ever be possible to 
simulate. Nevertheless, the results presented here show very 
little dependence of the critical value of C a on ReM- For 
Pi"m = 1, for example, we find C" 1 * = 1.2 for ReM ~ 6 
and C£ rit = 1.5 for Re M « 600. On the other hand, for 
larger values of Pi"m, the value of C£ rit can drop below unity 
(C -crit = o g for PrM = 100 ) wijjjg these chan ge S f C^ rit 

are theoretically not yet understood, it seems clear that they 
are small and do not provide support for an entirely different 
scaling law, as anticipated in recent work ifllll . 
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